Fluctuating Topological Defects in 2D Liquids: Heterogeneous Motion and Noise 
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We measure the defect density as a function of time at different temperatures in simulations of a two 
dimensional system of interacting particles. Just above the solid to liquid transition temperature, 
the power spectrum of the defect fluctuations shows a 1// signature, which crosses over to a white 
noise signature at higher temperatures. When 1/f noise is present, the 5-7 defects predominately 
form string like structures, and the particle trajectories show a ID correlated motion that follows 
the defect strings. At higher temperatures this heterogeneous motion is lost. We demonstrate this 
heterogeneity both in systems interacting with a short ranged screened Coulomb interaction, as well 
as in systems with a long range logarithmic interaction between the particles. 
PACS numbers: 82.70.Dd, 52.27.Lw, 61.20.Ja, 61.72.Bb 
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In liquids and glassy systems there has been consid- 
erable interest in dynamical heterogeneities which occur 
when certain regions of the sample have a higher mo- 
bility than the rest of the sample. Particle motion in 
these systems often occurs by means of correlated mo- 
tion of a group of particles along a string like structure. 
0. Most of the studies of heterogeneities in 2D and 
3D have focused on systems that are supercooled near 
the glass transition in simulations || and experiments 
HHHI ■ A recent study considered a relatively simple sys- 
tem of a 2D dense monodisperse colloidal assembly which 
shows an ordered to disordered transition as a function 
of density [Q . In this system the colloids form a defect 
free triangular lattice at high densities and disorder for 
lower densities. Near the disordering transition at the 
close packing density, collective heterogeneous motion of 
the particles appears, and the system consists of domains 
of sixfold coordinated particles, n c — 6, surrounded by 
grain boundaries composed of strings of n c — 7 and 
n c = 5 dislocations. The formation of grain boundaries 
or defect condensation in 2D monodisperse liquids has 
also been observed in colloidal assemblies |5||(| and in 
dusty plasmas (fjjUIJ. 

Since the dislocations have a complex interaction with 
one another, consisting of a long range logarithmic re- 
pulsion and a short range attractive core interaction, one 
could expect complicated dislocation dynamics for sys- 
tems with even a low density of dislocations. In the liquid 
state, creation and annihilation of dislocation pairs also 
occurs due to the thermally induced motions of the un- 
derlying particles. Although there have been numerous 
studies of the fluctuations of the particles in 2D liquids, 
the fluctuations in the defect density, defined as the num- 
ber of particles with n c ^ 6, have not been investigated. 
This measure would be easy to access in experiments on 
systems such as colloids and dusty plasmas where the in- 
dividual particles can be directly imaged. It is not known 
how the density of n c ^ 6 particles would fluctuate as the 
ordering transition is approached. If the defects are con- 
centrated in clumps or grain boundaries, then it is likely 



that the creation or annihilation of defects will be highly 
correlated, which can give rise to 1/ f a fluctuation spec- 
tra. Conversely, if the defects are appearing and disap- 
pearing in an uncorrelated manner, a white noise spectra 
would arise. We also wish to connect the formation of 
strings of defects with the appearance of correlated par- 
ticle motion along ID strings. 

In this work we show that for 2D systems which form 
triangular lattices at low temperatures, for increasing 
temperature there is a transition from a non-defected 
regime to a defected regime where there is a proliferation 
of defects. Here we do not attempt to examine the na- 
ture of the disordering transition, such as whether there 
is an intermediate hexatic phase [lC[| . Instead, we con- 
centrate on the motions of the particles and the defect 
fluctuations in the liquid phase. Close to the disorder- 
ing transition, in the defected regime, the system consists 
of regions of particles with sixfold order surrounded by 
strings or grain boundaries of 5-7 fold defects, and there 
are few free dislocations. We also find that in this region, 
collective particle motion occurs in ID strings along the 
grain boundaries. When temperature is fixed, the den- 
sity and structure of the defect strings fluctuates with 
time. Near the disordering transition, the defect density 
shows large fluctuations with 1/f spectra. For increasing 
temperatures, the defect density increases and the string 
like structures of the defects and the ID string like mo- 
tion of the particles are lost. At the same time, the low 
frequency noise power increases and the spectrum crosses 
over to a white noise signature, indicating the loss of cor- 
relations in the creation or annihilation of defects. 

We numerically study 2D systems at finite tempera- 
ture using Langevin dynamics. We consider monodis- 
perse particles in a sample with periodic boundary con- 
ditions and have investigated two types of interactions. 
Most of the results presented here are for particles in- 
teracting via a Yukawa or screened Coulomb interaction 
potential, Vfcj) — (Q 2 /\ri — Vj |) exp(— n\ri — Yj\). Here 
r^) is the position of particle Q is the charge of the 
particles, 1/k is the screening length, and we use 
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FIG. 1. The Voronoi constructions for a 2D system of parti- 
cles interacting with a screened Coulomb potential. A particle 
is centered at each polygon and the polygon color is: n c = 6, 
white; n c — 5, dark gray, n c — 7, black, (a) T/T d = 1.04; (b) 
T/T d = 7.0. 



k = 2/ ao, where do is the lattice constant. We choose 
to study this interaction since the quantities we study 
in this work can be accessed in systems such as 2D col- 
loidal assemblies and dusty plasmas, in which the parti- 
cles interact via a screened Coulomb interaction. We also 
consider particles interacting via a logarithmic potential, 
V(r) — ln(r). Unlike the Yukawa potential, this inter- 
action is long-range so that an interaction cutoff cannot 
be implemented. For increased computational efficiency 
we use the summation method of Rcf. |ll[]. A physi- 
cal example of logarithmically interacting particles with 
overdamped dynamics is vortices in thin-film type-II su- 
perconductors. 

The equation of motion for particle i is: 

^ = _^V i F(r ij )+f T (1) 

where fx is a randomly fluctuating force due to ther- 
mal kicks with < f T {t) >= and < f T {t)i T (t') >= 



2kBTS(t—t'). In all the simulations we initially start from 
an ordered triangular configuration. We fix the temper- 
ature fx for 10 6 time steps before we begin to take data. 
We note that previous simulations have shown that ex- 
tremely long time transients can arise near the order to 
disorder transition of up to 2 x 10 6 time steps for sys- 
tem sizes of 36864 particles [Q. In this work we limit 
ourselves to system sizes of N c = 2000 or less, so that 
10 6 time steps is adequate for equilibration. In addition, 
histograms of the time series of the defect density are 
Gaussian, which is further evidence that our systems are 
equilibrated. We measure temperature in units of T/Td 
where is the temperature at which the first free dislo- 
cations appear. 

In Fig. 1(a) we show the Voronoi construction of a sys- 
tem of particles with a screened Coulomb interaction for 
T/Td = 1-04. The Voronoi construction is similar to the 
Wigner-Seitz cell construction and an individual particle 
is located in the center of each cell. If a particle has six 
nearest neighbors, n c — 6, then the polygon has six sides. 
In Fig. 1 particles with n c — 6 are white, n c — 5 are dark 
gray, and n c = 5 are black. We analyze a series of such 
images at fixed T/T d = 1.04. We find that 29% of the 
particles are defected and that most (94.5%) of the de- 
fects are part of a cluster or condensate rather than free. 
In some regions, the defects form strings or grain bound- 
aries. There are also a small number of free disclinations 
present, indicating that we are in the liquid phase rather 
than a hexatic phase, according to 2D continuous melt- 
ing theories. Clustering of defects has been observed in 
experiments on dusty plasmas which find compara- 
ble numbers of free dislocations. In Fig. 1(b) we show 
the same system at a higher temperature, T/Td = 9. 
Here, the number of defects is higher with most of the 
defects again predominately appearing in clusters; how- 
ever, there is no string like characteristic to the clusters. 
Instead, the defects form clump like objects. 

In Fig. 2(a) we show the positions of the particles 
(black dots) and trajectories (lines) for the same system 
in Fig. 1(a). The trajectories are taken for the particle 
motion over 10000 time steps. Fig. 2(a) shows that the 
motion is heterogeneous with certain regions in which 
the particles do not move. These areas are also the de- 
fect free regions. In the areas where there is motion, the 
particles move by a, and there is some tendency for the 
motion to occur along string like paths. We also often 
find places where the motion occurs in a circular path 
with an immobile particle in the center. A similar kind 
of motion in ID like strings and rotations was seen in 
experiments on dense colloidal liquids H . If we take tra- 
jectories for longer times, eventually all the particles take 
part in some motion. In Fig. 2(b) we show the positions 
and the trajectories for the high temperature phase in 
system Fig. 1(b) for 2000 time steps. Here the motion 
occurs everywhere with no evidence for correlated sta- 
tionary regions. We find that even for very short times, 
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FIG. 2. Particle positions (black dots) and trajectories 
(black lines) for a fixed period of time for (a) T/T d — 1.04 
and (b) T/T d = 7.0. 

there are no correlations in the motion. 

We next consider a way to characterize the behaviors 
of the system at the different temperatures by measur- 
ing the fluctuations of the defect density as a function of 
time. We compute the defect configuration and density 
every 20 time steps for 20000 frames and obtain a time 
series of the defect density vs time for several temper- 
atures from T/Td = 1-04 to 10. We have considered a 
variety of sampling rates and find consistent results. In 
Fig. 3(a) we show a portion of the time series of the frac- 
tion of n c — 6 particles, Pe(t), for T/Td = 1-04 (upper 
curve) and T/T d = 7.0 (lower curve). For T/T d = 1.04 
the fluctuations show long time variations, while for the 
higher temperature the fluctuations are very rapid. In 
Fig. 3(b) we plot the power spectrum S(f) of Pe(t) for 
T/Td = 1-04, which fits well to a 1/ f a scaling over more 
than three decades with the best fit a = 1.04, close to 
1/f noise. As the temperature increases, the spectrum 
changes from 1/f to white noise (a = 0) for low frequen- 
cies. The small frequencies correspond to the long 
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FIG. 3. (a) A portion of the time series of the fraction of 
n c — 6 particles, Pe(t), for (upper curve) T/Td = 1-04 and 
(lower curve) T/T d — 7.0. (b) The power spectrum for the 
time series at T/T d = 1.04. The solid line indicates a slope of 
1/f. (c) The power spectrum for the time series at T/T d = 7.0 
along with a 1// line. 

time correlations in the system indicating that correlation 
times of the defect creation or annihilation decrease at 
higher temperatures. In Fig. 3(c), the power spectrum 
for T/Td = 7.0 is white with a = for a large portion of 
the curve. We have also considered Pe(t) for T/T d < 1.0 
where there are a small number of bound dislocation pairs 
present. Here we find white noise with small amplitude 
fluctuations. 

The magnitude of the noise power, obtained by inte- 
grating S(f) over the first octave of frequencies, is related 
to the defect density and the disordering transition. In 
Fig. 4(a) we show P e as a function of T and in Fig. 4(b) 
we show the corresponding noise power So. The peak 
in So coincides with the onset of the defect prolifera- 
tion. For increasing T, the noise power falls and saturates 
when the spectrum becomes white. There is a finite So 
for T/Td < 1-0 since, as noted, pairs of bound disloca- 
tions can still be thermally created. We have repeated 
the same set of simulations for particles interacting with 
a long range logarithmic potential and again find a noise 
power peak coinciding with the defect proliferation 



3 




2x10' 



T/T\ 



FIG. 4. (a) The fraction of six-fold particles Pe, vs T for 
a system of particles with a screened Coulomb interaction. 
Inset: Pq vs T for a system with logarithmic interactions, (b) 
The integrated noise power So, obtained from S(f), vs T for 
the same system in (a). Inset: So vs T for a system with 
logarithmic interactions. 

(insets of Fig. 4) and a 1// power spectrum just above 
T d . 

We have also considered a sudden quench from T/Td> 
1.0 to T /Td < 1.0. Here the system is out of equilibrium, 
and approaches equilibrium by means of the annihilation 
of the defects. In Fig. 5 we show the particle motions 
after a quench from T/T d = 1.04 to T/T d = 0.2. We 
find that the number of defects decreases rapidly to a 
saturation point where a few defects or grain boundaries 
remain. The motion of the particles corresponds to the 
annihilation of the defects. As seen in Fig. 5, the motion 
during the defect annihilation process has the same string 
like nature as the equilibrium T > T d motion near the 
defect proliferation transition. This result suggests that 
the heterogeneous particle motion is produced by the mo- 
tion of defects, particularly the creation or annihilation 
of these defects in a correlated manner. 

In summary, we have proposed a new way to exam- 
ine dynamical heterogeneities in a liquid by measuring 
the fluctuations in the topological defect density. We 
find that near the onset of defect proliferation, the par- 
ticle motion is heterogeneous and the defects cluster to- 
gether into grain boundaries or strings. The defect fluc- 
tuations in this phase have a 1// character and a large 
noise power. For higher temperatures the heterogeneities 
are lost and the fluctuation spectrum is white. We have 
also examined the defect annihilation after a quench from 
the disordered phase to the ordered phase and find that 
the particle motion also shows heterogeneities as the de- 
fects are annihilated. We have considered systems with 
screened Coulomb interactions and logarithmic interac- 
tions and find similar behavior. Our predictions can be 
easily tested in 2D dense colloidal liquids. It would also 
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FIG. 5. The particle positions (black dots) and trajecto- 
ries (black lines) for a system quenched from T/Td = 1.04 to 
T/Td = 0.3. The trajectories are analyzed once the system is 
at the lower temperature. During this time a portion of the 
defects are annihilated. 

be interesting to study the fluctuations of defects in 3D. 
This would be experimentally possible in 3D colloidal 
assemblies using confocal microscopy. 
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